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Abstract 

We investigate the dynamics of a scalar field governed by the Lifshitz-type action which should appear 
naturally in Horava-Lifshitz gravity. The wave of the scalar field may propagate with any speed without an 
upper bound. To preserve the causality, the action cannot have a generic form. Due to the superluminal 
propagation, a formation of a singularity may cause the breakdown of the predictability of the theory. To 
check whether such a catastrophe could occur in Horava-Lifshitz gravity, we investigate the dynamics of a 
dust. It turns out that the dust does not collapse completely to form a singularity in a generic situation, 
but expands again after it attains a maximum energy density. 



^e-mail address : suyama@phya.snu.ac.kr 



1 Introduction 



The construction of a quantum theory of gravity is one of the long-standing problems in theoretical 
physics. The main difficulty is the non-renormalizability which can be easily seen from the fact that the 
Newton constant has a negative mass dimension. Recent investigations (see [T] and references therein) on 
the perturbative amplitudes in supergravity theories seem to indicate that there would exist an unknown 
mechanism of cancellation among diagrams which might make some supergravity theories renormalizable. 
There is another proposal, called Horava-Lifshitz gravity [2] [2], which has attracted much attention 
recently. The basic idea of the approach is to modify the UV behavior of general relativity so that the 
perturbative renormalizability can be realized. The price to pay for this achievement is the abandonment 
of the Lorentz symmetry in the high energy regime: in this context, the Lorentz symmetry is regarded 
as an approximate symmetry observed only at low energy. See also [36j [37] [38] [39] [40] [41] . 

Since the Lorentz symmetry does not exist in the high energy regime, the notion of causality should 
be largely modified. For example, since the dispersion relation of the propagating modes is modified, 
they can propagate with any speed. This seems to raise a puzzle if this superluminal propagation is 
realized in the presence of a singularity. It would be natural to expect that a singularity would produce 
particles with an arbitrarily high energy. If they are produced, they can propagate with arbitrarily 
high speed, and therefore, there would always exist a particle which can reach an observer at infinity. 
Since such a particle would carry information of the singularity, the possibility of this phenomenon 
shows that the predictability of the theory breaks down. In general relativity, it is believed that any 
singularity should be hidden behind the event horizon, and it cannot affect any physical phenomena 
outside the event horizon. However, it does not seem to be the case in Horava-Lifshitz gravity since the 
light cone structure is not relevant, and therefore the "event horizon" defined as in general relativity 
cannot hide the singularity. The black hole solutions in Horava-Lifshitz gravity were investigated in 
[4][5][6][7][S][9][I0][II][I2][13][I4][15][16][I3 The other subtle issues in Horava-Lifshitz 

gravity were pointed out in [31] [32] [33] [34] . The dynamics of an anisotropic classical mechanical system 
was discussed recently [44] . 

In this paper, we investigate the propagation of a scalar field whose action is of the Lifshitz-type. We 
employ the WKB approximation, and study the trajectory of the ray of the wave of the scalar field. As 
is expected, it turns out that the wave can propagate with an arbitrary speed. Because of this, the light 
cone structure cannot have a physical relevance. Therefore, in Horava-Lifshitz gravity coupled to matter, 
there might occur the breakdown of the predictability, as pointed out above. To check this, we investigate 
a gravitational collapse in Horava-Lifshitz gravity, and check whether a singularity would be formed. We 
choose the simplest situation: the collapse of a dust preserving the homogeneity and the isotropy. Naively, 
this would be the most probable situation for a singularity to fornj§. Interestingly enough, it turns out 
that the dust does not collapse completely in a generic situation. Instead, at a maximum energy density, 
the dust ceases to collapse, and then starts expanding. The repulsive force necessary for this behavior is 

In this analysis, we assume, just for simplicity, that the equation of state is always the same. It may be expected 1431 
that the equation of state may change with the energy scale, but the inclusion of this property to the analysis will not be 
compatible with the homogeneity and the isotropy, and therefore, the analysis will become complicated. It will be very 
interesting to study this issue further. 
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induced by the higher derivative terms in Horava-Lifshitz gravity which are necessary to modify the UV 
behavior. Note that the analysis may be closely related to [n]I12[lS]ES][lI][lHlES][3nj some of which 
discussed a bounce cosmology in Horava-Lifshitz gravity. It should be noted that the analysis in this 
paper cannot be a proof of the absence of singularities in Horava-Lifshitz theory. The investigation on 
situations more general than the dust collapse is beyond the scope of this paper. 

This paper is organized as follows. In section 1, we investigate the propagation of a scalar field 
governed by a Lifshitz-type action, using the WKB approximation. Section 2 argues the issues related 
to causality. A gravitational collapse in Horava-Lifshitz gravity is discussed in section 3. Section 4 is 
devoted to discussion. 



2 Geometric optics for Horava-Lifshitz scalar 



We consider a complex scalar field whose dynamics is governed by the action 



where 



(2.1) 
(2.2) 



Dt<j}:=dt(j}-N'd,<j). 

and i nms from 1 to 3. The fields N,N'^,gij are combined to form a four-dimensional metric 

ds^ = -N'^dt^ + Qijidx' + N''dt){dx' + N'^dt). (2.3) 

The expression f {(f) , (f>\W , g) is a scalar with respect to the three-dimensional diffeomorphisms. More 
explicitly, /(</>, V, (?) is assumed to have the form 



M 



■^■n]l■■■3■. 



(2.4) 



n=0 



where G'*^ ' '"^'^ '^" is a three-dimensional tensor constructed from gr*-' , and is the covariant derivative 
with respect to gij. In the following, N is assumed to be independent of the spatial coordinates . 
The equation of motion of cf) is 



1 



dt 



M 

— V.(7V^A0)-5](-l)"G''^ 



■«r.Jl---J„V7. ... V- V ■ 



Vj„^ = 0. (2.5) 



It may be possible to analyze the propagation of 4> directly. However, the analysis will be easier if we 
restrict ourselves to deal with "geometric optics" , that is, a particle mechanical description of the ray of 
(j>. We use WKB approximation to obtain such a classical mechanical system. Assuming — e*'^ with S 
real, and neglecting sub-leading terms, we obtain 

1 ^ 

— (A5)' - E G^---^--^-{d,,S) ■ ■ ■ (a,„5)(%5) • • • (9,„5) = 0. (2.6) 



n=0 
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The left-hand side can be regarded as a Hamiltonian H by replacing df^S with p^. Then the equations 
of motion can be obtained as the Hamilton's equations. For example, if we started with the relativistic 
action, then the resulting Hamiltonian would be 

i^rci = -5^%P.-m2, (2.7) 

and the Hamilton's equations with this Hamiltonian 

^ = -23^"-^., ^ = d,g'"^p,p., (2.8) 

are equivalent to the geodesic equation with respect to the metric g^^. Note that the geodesic equation 
can also be derived from the conservation law of the energy-momentum tensor. 
For the general case, the Hamiltonian is 

H^^{pt-N^P^?-h{p^), (2.9) 
where p^ — g^-'piPj and h{x) is a polynomial defined as 

h{p^)^ f {1,1, p,g). (2.10) 
For the small momentum region, we expect that the system is approximately relativistic. This implies 

/i(p2) ^m2+p2+o(/). (2.11) 

We require that the (Euclidean) path-integral for this system should be well-defined. This implies 



lim h{p^) = +00. (2.12) 

— >oo 



The Hamilton's equations are 



r/r^ 2 

— = - — ip,-N^p^)N^-2h'{p')g'^p,, (2.13) 

dt 2 , , ^ , , 

- = —ip,-N%), (2.14) 



dT 



[pt ~ N^Pj)pkd,N'' + h'{p')d^g'''pkPu (2.15) 



^ - ^{pt-N^Pi)PkdtN^ + ^dtN{pt-N^p,f + h'{p')dtg^'pm- (2.16) 



This shows that a ray of does not follow any geodesic defined by the metric (|2.3p . 

In general, it is very difficult to eliminate the momenta using the velocities in the above equations. 
Let us consider a simple case in which the h{x) is a monomial, that is, 

Hr. = ^{pt-N^p^f-^{p'r (2.17) 
with a positive c. In this case, we obtain 

v' = -^{Pt^N^P^), (2.18) 

v' + N'v' = -cip^T'^g'^Pj. (2.19) 
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The momenta are given explicitly in terms of the velocities, and we obtain the Lagrangian 

L = v*pt + v'p, - H 

= —{v^f-c'iV') — , (2.20) 

where 

c'=(l-i-)c-^, (2.21) 

and 

y2 = g.jV'V^, V ^v' + N'v*. (2.22) 



3 Causality 



Since the Hamiltonian (|2.9p is different from the relativistic one (|2.7p . the equations of motion are 
different from the geodesic equation for the metric (|2.3p . Therefore, the issue on the causality would be 
very different from the familiar relativistic one. 

In our setup, the existence of the absolute time is assumed. This implies that the presence of closed 
time-like trajectories is forbidden. In this case, a typical causality-violating trajectory is the one which at 
first moves forward in time but at a later time moves backward in time. In other words, there could be a 
closed trajectory. Such a trajectory is possible only if ^ = is realizable. This implies that pt — N^Pi = 0, 
and therefore h{p'^) = must be realized at some value of r. As a result, such a causality violation never 
occurs if h{p'^) is always positive. (This is of course the case for the relativistic system.) On the other 
hand, if h{p^) allows to have a zero, then, in general, there would exist a trajectory which "turns around" 
at some value of t since is non-zero in general at the time when the trajectory turns around. 

Let us examine the simple case in which = 1 and A^' = 0. In this case, a trajectory turns around 
when pt = Q, and at this moment 

^^2^^ 2h'{j>^)d,9^Wpi- (3.1) 

The value of the right-hand side depends on the background metric, and is generically non-zero unless 
the background is static. 

There is another possibility of a strange behavior which never happens in the relativistic system. 
One notices that the expressions in the right-hand side of the Hamilton's equations are simplified when 
h'{p'^) = is realized. At this point in the phase space, the correspondence between the momenta and 
the velocities is lost. For example, consider again the case A^ = 1 and A^* = 0. Then, the r-derivatives of 
a;'. Pi and pt all vanish when h'{p^) — 0, and ^ — 2pt. In addition, one can show that 

= 2ptdtg'''pwi- (3.2) 

Therefore, if the background is static, then p^ and also pt are constant, due to the conservation of the 
Hamiltonian, and 

= const. X* — const. t — 2ptT + to (3.3) 
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is a solution. It looks very strange since non-zero pi correspond to the trajectory at rest as long as 
h'{p'^) ~ 0. Note that in a generic background, since -^p^ would not be zero when h'{p^) — 0, and 
therefore, such a strange solution is not allowed. 

In the following, let us assume that h(x) is a monotonically increasing function with h'(x) > in 
order to avoid the strange behaviors of the trajectory mentioned above. In this case, the causality is 
kept in the sense that there is no trajectory which moves backward in time, and therefore, no closed 
trajectory. However, the notion of the causality is still quite different from the one in general relativity. 
In particular, as is easily anticipated from the action (|2.ip with which we started, there is no physical 
relevance of the light cone. In fact, one can calculate the relativistic norm of the velocity as follows: 

g^^vf'v'' = -N^-{vy+gij{v' + N'v'){v^+N^v') 

- A[^h{p^)+p^{h'{p')n (3.4) 

This norm is constant only when the action (|2.ip is relativistic. This becomes indefinitely large and 
positive when the spatial momenta pi become large. This implies that a trajectory which was time-like at 
first may becomes space-like, and therefore, the trajectory can escape from the light cone defined by the 
metric (|2.3|) in the usual manner. Note that this property of the trajectories is inevitable as long as the 
degree of h{x) is larger than 1, that is, as long as the UV behavior is modified. In contrast, the presence 
of a closed trajectory, for example, can be avoided by choosing a suitable functional form of h{x). 

The irrelevance of the light cone poses the following question: is the formation of a black hole in 
Horava-Lifshitz gravity allowed? Suppose that black holes are defined in the same manner as in general 
relativity. In Horava-Lifshitz gravity, the event horizon cannot hide the singularity inside, since there 
exist trajectories which can escape from the inside of the event horizon. Therefore, it seems that the 
predictability of the theory would break down whenever a singularity is formed. In the next section, we 
will investigate the collapse of a dust in Horava-Lifshitz gravity. We choose this situation since, due to 
the absence of the pressure, the formation of a singularity by the collapse of the dust seems to be easily 
realized. Interestingly enough, it will be shown that the higher derivative terms prevent the dust from 
collapsing completely, and the dust will finally form a compact object with a finite size. 



4 Dynamics of dust in Horava-Lifshitz gravity 



Let us first recall the analysis of the dust collapse in general relativity [23] . Since Horava-Lifshitz gravity 
typically includes a cosmological constant, we include it to the analysis in [23] so that the effects of the 
higher derivative terms may be highlighted. 
We consider the following action 



S 



I 



d X\J —g 



-^R + a + Lr, 



(4.1) 



where = 327rG. For simplicity, we consider the collapse which preserves the homogeneity and the 
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isotropy. Then the metric is of the Robertson- Walker form 



-dr + a{tyhij{x)dx^dx^ , 

KAx)dx'dx^ = — + r^ide^ + sin^ Odcf)'^). 

1 — kr-^ 



(4.2) 



Let p{t) be the energy density of the dust. The conservation law determines p — poa ^. The Einstein 
equation reduces to 

'dt ) ~ 12 



Po 2 
era 

a 



(4.3) 



where we rescaled the radial coordinate r so that a(0) = 1 is satisfied. We assume that the collapse starts 
at i = 0, that is, we choose the initial condition d(0) = 0, which fixes 

k^^ipo~a). (4.4) 

In the case a — 0, there is an exact solution of (j4.3p . The solution is a cycloid whose parametric expression 
is 



i(l + COSV'), 
1 

{'ip + sin V')- 



2y/k 



The collapse a = 0, at which p diverges, occurs at a finite time t = T where 



T = - 



kVp(O)' 



(4.5) 
(4.6) 

(4.7) 



The behavior of a near a = is similar to the above exact solution even for a nonzero cr since its 
contribution is negligible when a is small. 

Next, we consider the same situation in the context of Horava-Lifshitz gravity. Following [24], we 
consider a rather general action 



S 



d^xdty/gN 



-CRijR'^ + + ^R + a + L„ 



where Lm is the matter Lagrangian. Various quantities are defined as 

K^J = ^(.gy-v,7V,-v,7V,), 

K = K,,g'^, 

Akl / 1 

^/9 \ 4 



(4.8) 

(4.9) 
(4.10) 
(4.11) 



and Rij , R are the Ricci tensor and the scalar curvature with respect to g^j . 

As in the relativistic case, let us consider a dynamics of the dust which preserves the homogeneity 
and the isotropy. We choose the Robertson- Walker metric 



ds^ = -N{tfdt^ + a{tfh,j{x)dx^dx^. 



(4.12) 
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Here we consider the so-called projectable theory. The spatial metric gij = a?hij satisfies 



2ka ^gij. 



(4.13) 



Therefore, the covariant derivative of the Ricci tensor vanishes. In addition, since hij is conformaUy 
flat, the Cotton tensor C*-' for hij vanishes. Inserting the metric ansatz into the action, we obtain the 
foUowing reduced action: 



Sr = V dt^N 



]^(3-9A)(- 



12fc2(C + 3?7) 6fc^ 



(4.14) 



We assume that the matter consists of a dust. The coupling of the dust to the metric is determined 



by 



where 



1 

2 , 

i / d-'xdt^N SN{-2NT'"') + 5N\2g,jN3T'' 



The equation of motion of N is 



3a(l-3A) - 



12fc2(C + 37?) 6fcC 
7a ^ 



a - p. 



(4.15) 
(4.16) 

(4.17) 



Here we chose iV = 1. 

A remark on the derivation of (|4.17p is in order. Since we assumed that N depends only on t, the 
equation of motion of N does not provide a local equation in general. Instead, one obtains 







(4.18) 



In the present case, however, the spatial dependence of the integrand disappears due to the ansatz (|4.12|) . 
and therefore, the integral provides just a trivial volume factor. In this way, one can obtain an equation 
which looks local, enabling one to obtain the above Friedmann-like equation (|4.17p . 

The matter density p is given in terms of a through the conservation law. The invariance under the 
time reparametrization provides 

(4.19) 

Note that, as in the case for the equation of motion of N , the conservation law for the time reparametriza- 
tion is not a local equation in general, since Horava-Lifshitz gravity preserves only a restricted part of the 
diffcomorphisms. In the case here, the ansatz of the homogeneity and the isotropy enables one to obtain 
(|4T9l) . The equation (jiTQ]) implies 



Pit) 



(4.20) 
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As in the relativistic case, we rescale the radial coordinate r so that a(0) — 1. 
Now the equation (|4.17p becomes 



-F(a), (4.21) 



da 

Via) = -^-^-co-C2a2, (4.22) 
a 



where 



4fc^(C + 3r/) 

= a(l-3A) ' ^^•2^) 

- -RT^A)' (^-^^^ 

We consider a generic case where ( + 31] ^ 0. This is the case when the theory satisfies the detailed 
balance condition, and a small deviation from the detailed balance may not change the sign of C + 3ry. 
We choose the initial condition d(0) = 0, a(0) = 1 as in the relativistic case. This implies 

C-2 + C-1 + Co + C2 = 0, (4.27) 

which determines k. This equation is second order in k, while k is assumed to be real. A real k is allowed 
iff 

4(C + 3r?)(po-a) + 3e' >0. (4.28) 

In general, there is a singular situatiorl§. If po = is satisfied, then fc = is a solution of (|4.27p . In 
this case, c_2 and Cq in V{a) vanish. As long as A > i, the resulting equation is the same as (14. 3p up 
to a rescaling of t, implying that the dust collapses completely and a singularity is formed. Note that 
the formation of the singularity in this case should be exceptional due to the highly symmetric setup we 
assumed. The solution with fc = has the spatial hypersurface which is flat, and the higher derivative 
terms do not contribute to the dynamics of a. This is due to the cancellation between the positive energy 
Po and the negative energy —a. In general, a small deviation from the flat spatial hypersurface would 
induce a small but non-zero c_2, at least effectively, and therefore, the behavior of the solution would be 
similar to the one with fc ^ which will be investigated below. 

Let us consider the case in which the detailed balance condition is satisfled. In this case, the coefficients 
in the action (|4.8p are parametrized as follows [M] : 

a = 4, (4.29) 

C = (4.30) 

1 - 4A 

^ ^ 8(1 -3A) 4 ' (^-'^^ 



We would like to thank Soo-Jong Rey for pointing out this singular situation. 
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Figure 1: The plot of v{a) = ^^^^^V{a) with fc = i. 



K-2//2 



.2 2 2 

They give, for example, c_2 = — xg(3A-i)^ • Using those parametrization and after a suitable rescaling of 
t, the equation (|4.2ip becomes 

fdaV (k V 8(3A-l)pol 

where k — The condition (|4.28p for the existence of a solution reduces to A > i, and the solutions 
of (14281) are 

In the following, we concentrate on the case A > i. Typical shape of V[a) is depicted in Fig. [TJ 

It turns out that V{a) is a concave function of a, and the equation V{a) = has two real solutions. 
We expect that the dust would contract at first, and this is compatible with (|4.34p only if F'(l) > 0. 
The value of y' (1) is 

^'(l)-Y^^f^(fc + 3)(l-fc), (4.36) 
and this quantity can be positive iff — 3 < fc < 1, implying 

PO < Pmax = 3X^1 ■ ^'^■^'^ 

This condition shows that there exists a maximum possible energy density determined by the parameters 
of the theory. If the initial density is above the maximum, the dust would start expanding rather than 
contracting. 
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0.6 



0.2 



Figure 2: The time evolution of the scale factor a{t) for k — ^. 

The qualitative behavior of the solution of (|4.2ip can be easily understood through the analogy with 
the corresponding classical mechanical system in which a is a coordinate variable: 



—j +Via)^0. (4.38) 

Suppose that po < pmax is satisfied. The value of a starts decreasing at first, due to the gravitational 
interaction, and the energy density becomes large. However, the term ^ in (|4.34p . which is absent in 
the relativistic case (I4.3p . acts as a repulsive force, and it prevents the dust from collapsing completely. 
Instead, the system will attain a minimum value of a at some finite time, and then expand again. The 
expansion will also stop in a finite time, and the system will oscillate with a finite period. A numerical 
solution for fc = i is shown in Fig[2l 

Note that the appearance of this periodic behavior would be due to the fact that we have considered 
a highly symmetric situation. If there would be a small asymmetry, then the oscillation of the dust 
would probably emit the gravitational wave, and the system would settle down to a static configuration; 
a compact ball of the dust with a finite size. The energy density would be also finite during the process. 

It is fairly obvious that the qualitative behavior of the dust discussed above is preserved by a small 
perturbation of the theory away from the detailed balance, as long as c_2 is kept negative. For example, 
the addition of a term proportional to R simply shifts the "energy" of the mechanical system (I4.38P , which 
will not change the behavior of the solution as long as the condition V{a) = has two real solutions for 
a. Another simple deviation from the detail balance case can be achieved by parametrizing a as 

with another parameter a. The equation (j4.34p becomes 



daV (~k \ 8(3A-l)pol 



dt J \a j 3 K?iJ?h? a 



da\ (4.40) 
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Figure 3: The plot of v{a) with a — I and k — ^. 



It turns out that the behavior of the solution is qualitatively the same as the one for ct = 0, as long as 
(7 < 1 is satisfied. Let us consider the case ct = 1 in which the cosmological constant cr vanishes. The 
initial condition for a determines k as 



which is compatible with the condition V'{1) > 1, or equivalently, —2 < A; < 0. The shape of V{a) 
with A; = — 1 is depicted in Fig. [31 It is concluded that even in the case a = 1, there is a solution 
for the dust which is oscillating. Notice that there is a crucial difference from the result in general 
relativity, in addition to the absence of the singularity formation. The value of k is always positive in 
general relativity without cosmological constant, but in Horava-Lifshitz gravity k is always negative if 
the cosmological constant vanishes. This difference may affect the space-time geometry outside the dust 
through the continuity of the metric at the surface of the dust. Therefore, the difference between general 
relativity and Horava-Lifshitz gravity might be observed by an observer at infinity. If a is sufficiently 
large, there is no solution in which the dust starts contracting at first. This is because of the expansion 
due to the large positive cosmological constant. 

We have considered the dynamics of the dust satisfying the initial condition a(0) = l,a(0) = 0. It 
is very easy to generalize the above analysis to the case in which a(0) = is still assumed but a(0) can 
be arbitrary. This generalization is necessary if one would like to consider a realistic situation where the 
higher-derivative terms are small at the beginning of the collapse. 

Let us define the normalized scale factor a. 




(4.41) 



a{t) : 



ait) 



(4.42) 



a(0)- 
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In terms of a, the equation (|4.40p becomes 

fd-aV ( ~k/a{Or X , 8 (3A-l)po/a(0)3 ^ 

The effect of the generic initial value of a{t) appears as a simple rescaling of k and po. It is easy to see that 
the qualitative behavior of the dust does not change by varying the value of a(0). If a(0) would be taken 
to be large, then fc/a(0)^ would be neghgible at first, and the behavior of the dust is well approximated 
by the collapse in Einstein gravity with a negative cosmological constant reviewed at the beginning of 
this section. When a becomes small enough so that the effect of k becomes relevant, the dynamics of the 
dust deviates from the complete collapse in Einstein gravity, and turns into expansion. 



5 Discussion 



We have discussed the dynamics of a matter in Horava-Lifshitz gravity. It would be natural to expect that 
the modification of the UV behavior of the gravity part should be accompanied by a similar modification 
for the matter part. We investigated the propagation of a scalar field governed by such a modified action 
in the WKB approximation. The consistency of the scalar propagation with the causality restricts the 
form of the action. For the consistent scalar field, there is no relevance of the light cone. Therefore, 
in contrast to general relativity, a formation of a singularity would directly imply the breakdown of the 
predictability. We also investigated the dust collapse in Horava-Lifshitz gravity. We found that any 
singularity would not be formed during the dust collapse preserving the homogeneity and the isotropy 
which is quite different from the situation in general relativity. It is tempting to speculate that the 
formation of any singularity is forbidden in Horava-Lifshitz theory, but the investigation on this issue is 
left for further research. The breakdown of the predictability due to a singularity formation would not 
be realized in Horava-Lifshitz gravity if the formation of the singularity would be forbidden. 

The repulsive interaction which prevents the dust collapse is induced by the higher derivative terms. 
However, it is clear that every higher derivative terms do not necessarily provide such a repulsive inter- 
action. In fact, one can choose the action so that c_2 in (|4.2ip is positive. The negative sign of c_2 
should be related to the fact that the "potential term" of Horava-Lifshitz gravity is positive definite in 
the detailed balance case. Due to this structure, a high-curvature configuration costs much energy than 
a low-curvature one, and therefore the former is not favored. It seems that the choice of c_2 might 
be restricted since c_2 could be related to the high energy behavior of the graviton propagator. It is 
interesting to investigate this point further. 

One interesting issue to be studied next would be the outer solution for the dust collapse. In this 
paper, we only considered the solution which describes the inside of the dust. If the dust is a ball with a 
finite radius, then the metric describing the outside of the ball would be spherically symmetric but is not 
homogeneous. In general relativity, one can find the solution describing the outside of the dust, thanks 
to Birkhoff theorem. It is very interesting to study whether there is such a theorem in Horava-Lifshitz 
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gravity. Since the available diffeomorphisms in Horava-Lifshitz gravity are less than those in general 
relativity, such a theorem, if exists, must be a much weaker one. 

Another interesting issue would be the gravitational radiation in Horava-Lifshitz gravity. As we 
mentioned at the end of the previous section, the construction of a more realistic model of the dust 
collapse might require the understanding of the emission of the gravitational wave from the oscillating 
matter. From the physical ground, the spherically symmetric wave should be forbidden also in Horava- 
Lifshitz gravity as long as it is assumed to approach general relativity in low energy. Since otherwise the 
emitted gravitational wave at a high energy region near the center may propagate toward the outer low 
energy region where such a wave cannot exist as in general relativity. If a spherically symmetric wave 
does not exist, then the variety of the spherically symmetric solutions would be rather limited, and a 
weaker version of Birkhoff theorem could be expected. It is interesting to clarify this issue. 
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